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SUM^RY 


The literature contains several analytical and experimental evalua- 
tions of Mode I stress-intensity factore for corner cracks at holes in 
plates subjected to remote tension, remote bending, or pin loading in 
the hole. Unfortunately, these solutions give very different stress- 
intensity factors for the same crack configuration and loading. 

The pui*po8e of this paper is to present stress-intensity factors, 
calculated by a three-dimensional finite-element analysis, for shallow 
or deep quarter-elliptical corner cracks at the edge of a hole in a finite- 
thickness plate. TYie plate was subjected to remote uniform tension, 
remote bending, or simulated pin loading in the hole. A wide range of 
configuration parameters wan investigated. The crack depth-to-plate 
thickness ranged from 0.2 to 0.8, while the ratio of crack depth to crack 
length ranged from 0.2 to 2. The ratio of hole radius to plate thickness 
was held at 0.5. To verify the accuracy of the three-dimensional finite- 
element models employed, convergence was studied by varying the numbers 
of degrees of freedom (the number ranged from UUOO to 9300). The stress- 
intensity factor variations along the crack front are compai*ed with 
solutions from the literature. 
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Corner cracks at holes are among the most common flavs in aircraft ' 
structures. Accurate stress analyses of corner-crack configurations are 
needed for reliable prediction of crack-growth rates and fracture strengths, 
However, because of the complexities of such problems, exact solutions are 
not available. All investigators have used experimental or approximate 
analytical methods to obtain stress-intensity factors. 

Engineering estimates for Mode I stress-intensity factors for some 
corner-crack configurations have been made by Hall and Finger [l], Liu [2], 
and Newman [3]. These investigators did not consider the variation of 
stress-intensity factor along the crack front. Rather, their estimates 
gave a single value of stress-intensity factor for each crack configuration 
and, therefore, their estimates might be considered only as an average * 
value of stress-intensity factor along the crack front. Shah [H] used the 
alternating method, along with an engineering estimate, to calculate the 
stress-intensity factor variation along the crack front for the corner- 
crack configuration subjected to either remote tension or pin loading in 


the hole. 
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McGowan and Smith [5], Sjoith, Jolles and Peters [6], and Smith, 

Peters and Gou [7]» used three-dimensional photoelastic techniques to 
obtain stress-intensity factors for a variety of eorner-craclc configu- 
rations subjected to either remote tension or pin loading in the hole. 

They reported stress-intensitjr factors at two or three locations along 
the crack front. 

A few three-dimensional stress analyses of the corner' crack con- 
figuration have been reported recently. To analyze a quarter-elliptical 
corner crack emanating from a hole in a finite-thickness plate, Ganong [8] 

j! 

used the alternating method to obtain the- stress -intensity-factor varia- 
tions along the crack front, while Katbiresan [9l and Heckmer and Bloom [10] 
used three-dimensional finite-element methiods. Unfortunately, some of 
these solutions for the same configuration and loading gave widely varying 
stress-intensity factors. 

The purpose of this paper is to present Mode I stress-intensity 
factors, calculated by a three-dimensional finite-element analysis [11,12], 
for shallow or deep quarter-elliptical corner cracks at the edge of 
a hole in finite-thickness plates. The finite-thickness plates 
were subjected to remote uniform tension, remote bending, or simulated pin. 
loading in the hole. A wide range of configuration parameters was con- 
sidered. The ratio of crack depth to plate thickness ranged from 0.2 to 
0.8, th ratio of crack depth to crack length ranged from 0,2 to 2, while 

c\ 

the ratio of hole radius to plate thickness was held at 0.5. To study 
convergence , finihe-element models with to 9300 degrees of freedom 



were analyzed. The stresB~intensity factors were calculated "by using 
a nodal-force method [11]. The stress-intensity-factor variations along 
the crack front are presented and compared with other solutions from the 
literature. 


SYMBOLS 


a 

b 

c 


F. 


1 



Q 

R 


S 

P 


t 


x,y,z 



depth of corner crack 
half-width of cracked plate 
length of corner crack 

stress-intensity boundary-correction factor 
half-length of cracked plate 
stress-intensity factor (Mode l) 
applied bending moment 
applied load 

shape factor for an elliptical crack 
radius of hole 

remote outer fiber bending stress 
hole bearing stress (= P/2Rt) 
remote uniform tensile stress 
plate thickness 
Cartesian coordinates 
angular measuremei j (see Fig. 5 ) 


1 » 
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normal stress applied on hole boundary 
V Poisson's ratio 

<() parametric angle of the ellipse 

ANALYSIS 

The corner-crack configuration is shown in Figure 1. An elastic 
plate of thickness t, width 2b and length 2h, contains a through-the- 
thickness circular hole of radius R. Emanating from the hole are two 
symmetrically placed quarter-elliptical corner cracks of length c on 
the front surface and of depth a on the hole surface. 

Three types of loading have been considered: remote tension, remote 

bending, and wedge loading in the hole. The solution for wedge loading 
in the hole was also used In combination with the remote tension solution 
to obtain solutions for the case of simulated pin loading in the hole. 

ir 

I 

Finite-Element Idealization 

Two types of elements (isoparametric and singular [13]) were used 
in combination to model elastic bodies with;, quarter-elliptical corner 
cracks. Figure 2 shows a typical finite-element model for a corner crack 
in a large body, b/(R+c) > 5 and h/b = 1.8. This model idealizes one 
quarter of the body. Various n;xmbers of wedges were used to form the 
desired crack configuration in the y = 0 plane. A wedge is the region 
between two radial lines emanating from the point [x] = R, z = 0 in the 
y = 0 plane, for the circular crack model and is the region between two 


hyperbolas in the y = 0 plane for the elliptic crack model as shovm in v. 
Figure 3. For further details on modelling see Reference 12. The model 
shown in Figure 2 has eight wedges, each of which is composed of elements 
around the crack front that are identical in pattern to that shown in the 
x,y plane. The isoparametric (eight-noded hexahedron) elements were used 
everywhere except near the crack front. Around the crack front (such as in 
the x,y plane) each wedge contained eight ’’singularity” elements in the 
shape Of pentahedrons. The "singularity” elements produced a singular stress 
field at the crack front. Details of the formulation of these types of 

1 • 

elements are given in References 11 and 13 and are not repeated here. 

The finite -element model for the quarter-elliptical corner crack 
was obtained from the finite -element model for the quarter-circular 
crack by usiiig an elliptic transformation only in the region of the 
wedges. This transformation was convenient because the nodal forces 
along the normals to the crack front are used to calculate the stress- 
intensity factors. If (x,y,z) are the Cartesian coordinates of a node in 
the circular-crack model and (x',y',t’) are the coordinates of that same 
node in the elliptical -crack model, then the ti'ansformation is" 
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for X ^ R. Figure 3 shows how circular arcs and radial lines in the x,z 
plane of the circular-crack model are transformed hy equations (l) into 
ellipses and hyperbolas, respectively, in the x’ , a' plane of the 
elliptical-crack model. Because equations (l) are not valid at x = R 
and z = 0, a circle of very small radius, a/lOOO, was used near x = R 
and z = 0 in tbe circular crack model. The small circle, which maps onto 
a very narrow ellipse in the x' , z’ plane, avoids ill-shaped elements and 
allows the use of 8-noded elements in the elliptical-crack model. The 

elliptic transformation reduced the^b/c ratio; therefore, to eliminate the 

\ 1 

influence of plate width, additional rectangular prism elements were added 
along the x'-axis. 

Loading 

i> 

Three types of loading were applied to the finite -element models 
of the corner-crack configuration: (l) remote uniform tension, (2) remote 

bending, and (3) wedge loading in the hole. The loadings are illustrated 
in Figure U. The remote uniform tension is in Figure 4a; the remote 
outer fiber bending stress, S^, in Figure 4b is calculated from the applied 
bending moment; , and the applied normal stress on the hole boundary, cT^, 
in Figure 4c was assumed to be given by 

■Ji ? 

- C' !( 
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T 




where P is the total applied force acting in the y-direction over the arc 
from 9 = 0 to TT. The particular form of was chosen to simulate pin 
loading (see Fig. 5)* The inode I stress-intensity factor for the case of 
pin.J.oading in the hole is obtained by appropriate superposition of the 
results for remote uniform tension (Fig. 4a) and for wedge loading in'jhe 
hole (Fig. 4c). (See Ref. ' . 


Stress-Intensity Factor 

Only loadings which cause Mode I (tension) deformations were analyzed 
The Mode I stress-intensity factor, K^-, at any point along the quarter- 
elliptical corner crack in a finite-thickness plate was taken to be 





/, aT T, /a a R , i 

5 ^1 ‘t’ c- t’ 


( 3 ) 


where the subscript i denotes the type of applied loalliig (remote 
tension (i=t ) , remote bending (i=b), or wedge loading in the hole (i=p)); 
Q is the shape factor for an ellipse and is given by the square of the 
complete elliptic integral of the secondjkind [l5]- The boundary-correc- 
tion factor, F^, is a function of crack depth, crack length, hole radius, 
plate thickness, and the parametric angle of the ellipse. Note that the 
length and width of the plate were chosen large enougli so that they would 
have a negligible effect bn stress intensity. Values for Fj^ were 
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calculated as a function of ~ and (|) for ^ = 0.5* The ~ ratios 
ranged from 0.2 to 2 and the ^ ratios ranged from 0.2 to 0.8. 

The stress-intensity factor (Mode I) for the simulated pin-lo>d con- 
figuration shown in Figure 5 was obtained by superposition of the results 
for the configurations shown in Figures Im and tc. (Superposition of 

j; 

these results gives the correct normal stresses for the pin-load Configura- 
tion only along the y = 0 plane.) For this case, the Mode I stress- 
intensity factor is given by 



i. (S. p + p ) 
2 S t p^ 


(M 




where S = P/2Rt , and F. and F are the boundary -correction factors ob- 
P t p 

tained from equation 3 for remote tension and wedge loading in the hole, 
respectively. (This superposition method gives no information on Mode II 
stress-intensity __factorb. ) 

The stress -intensity factors from the finite-element models of "the 
quarter-elliptical corner cracks were obtained by using a nodal -force 
method, details of which are given in Reference 11. In this method, the 
nodal forces normal to the crack plane and ahead of the crack front are 
used to evaluate the stress-intensity_ factor . In contrast, the crack- 
opening displacement (COD) method ‘presented in [13] requires a prior 
assumption of either plane stress or plane strain, which is a potential 
source of inaccuracy. The nodal-force method requires no such assumption. 
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RESULT AMD DISCU’SSIOM 

U 

In the following sections , results for qyoi'ter~circialar and quarter- 
elliptical corner cracks emanating from a circular hole (R/t « 0.5) sub- 
jected to various loading are presented. Convergence of the stress- 
intensity factors for a deep quarter-elliptical corner crack was studied 
by varying the number of degrees of freedom in the finite-element models 
from to 9300. The stress -intensity factor variations along the crack 
front for quarter- circular (a/c 5= l) and quarter-elliptical (a/c = 0.2 

and 2) corner cracks are presented as functions of a/t with R/t = 0.5« The 

! 

stress-intensity factors are compared with other solutions and experimental 

if 

results from the literature. 

Convergence 

In reference 11 » the present method was applied to the problems of 
embedded circular (a/c - l) and embedded elliptical (a/c =0.2) cracks in 

a large bidy under uniform tension. Because the results obtained for 

''' // '' 

these crack shapes were generally within 1 percent of the exact solutions 
[15] » the present method is expected to be suitable for analyses of the 
more complex configurations considered herein, provided that enough degrees 
of freedom are used to ensure convergence. 

Figure 6 shows the results of the convergence study for a deep 
quarter-elliptical corner crack (a/c = 0.2, a/t = 0.8, R/t = 0.5 ). This 

if ■■ f , 
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fionriguratiori was chosen because proximity of the back surface to the 
crack front was expected to cause convergence difficulties. Four finite- 
element models with UHOO to 9300 degrees of freedom (DOF) were analyzed. 
These models used 2, H, 8 or 10 wedges to idealize the elliptic crack front 
on the y a 0 plane. A typical eight-wedge model is shown in Figure 2. 

The two-, four-, and eight-wedge models had equal wedge angles. The ten- 
wedge model had nonuhiform wedges with smaller wedges near the intersection 
of the crack with the hole surface {<{» « tt/ 2). The smaller wedges near the 
free surface were used to study the "boundary-layer" effect suggested by 
Hartranft .and Sih [l6]. The ideaJ.ization used to model the circular hole 
was the same for all models (see Fig. 2). The two finest models (81+54 and 
930o degrees of freedom) gave stress-intensity factors within less than one 
percent of each other except at the intersection of the crack front and the 
hole surface, where the difference was about 5 percent. The ten-wedge 
model with about 9300 degrees of freedom was used subsequently to obtain 
stress-intensity factors for all other crack configurations, 

Quarter-Circular Corner Cracks 

Figure 7 shows the normalized stress-intensity factors for a quarter- 
circular comer crack as a function of the parametric angle , (j> , and the 
crack-depth-to-plate-thickness ratio, a/t, for remote tension, remote 
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bending, or wedge loading In the hole. For all loading types considered, 
the smaller a/t ratios generally gave hi^er nonaalized stress«>lnten8lty 
factors. For renote tension and for wedge loading in the hole, the 
maximum normalized streBs^lntensity factor occurs near the intersecticxi 
of the crack with the hole surface ((^ » x/2). For remote bending, the 
normalized stress-intensity factor was largest at « 0 (front surface). 

3he negative stress-intensity factors shown in Figure 7b for a/t = 0.8 
are, of coarse, meaningful only in the presence of sufficient tensile 
loading to prevent contact between the crack surfaces. 

Ibe stress-intensity factors (Mode l) for quarter-circular comer 
cracks subjected to simulated pin loading, as shown in Figure can be 
obtained from the results shown in Figures 7a and 7c by use of equation b. 

Figures 0a and 8b show stress -intensity factors obtained by several 
investigators for a quarter-circular corner crack in a finite-thickness 
plate under remote tension. Figure 8a shows the results for a comer 
crack with a/t = 0.2 and R/t « 0.5. TSie present results are shown as 
solid circiilar symbols. Shah's estimate [U] is about 6 i>ercent hi^er 
than the present results. Ihe results from Heckmer and Bloom [10 ] are 
10 to 20 percent higher than the present results near the free surfaces 
(i^ * 0 and it/2), but only about 3 percent higher in the interior. Ibe 
photoelastic results [5,7lt which represent the experimental stress-intensity 
value at the approximate midpoint of each i^otoelastic slice, generally 
agree well with the present results. 
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The stress intensity is less near the intersection of the craclc with 
the hole surface {<{) = Tr/ti) than ^^ust below the sui'face. This behavior 
agrees qualitatively with that proposed by Hartranft and Sih [l6]. Note 
that the wedge angles near = Tr/2 are smaller than those near (j> = 0, 
Similar refinements in the wedge angles near <|) = 0 would probably have 
revealed lower stress-intensity factors in that region as well. 

Figure 8b shows the results for a deeper corner craJvt (a/t = 0.5) 
than that shown in Figure 8a. Shah's estimate [4] is within 10 percent 
of the present results, except that Shah does not predict the lower stress 
intensity near (j) = fr/2. Kathiresan's results [9] are 10 to 20 percent 
lower than the present results. This descrepancy may be attributed to the 
fact that the present study employed over three and one-half times as many 
degrees of fi*eedom as did Kathiresan's. Ganong's results [8], which are 
tVr a single corner crack emanating from a hole, are 25 percent higher than 
the present results. The results for a single corner crack are expected to 
be lower than those for two corner cracks [4], Again, the photoelastic 
results from Smith, Peters, and Gou [7] t^reed to within ±10 percent' with 
the present results. 

S-rrme of the discrepancies between the various finite-element solu- 
tions may be attributed to the method used to extract the stress-intensity 
factor from the three-dimensional analyses . Figure 9 shows the stress- 
intensity factor variation along the crack front for a typical corner-crack 
configuration as determined from the present analysis by two methods. The 
configuration and loading are those appropriate to Figure 8a. The stress- 
intensity factors were determined from the crack-opening displacements (COD) 
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asHuming either plane strees or plane strain and from the nodal forces ahead 
of the crack front [ll]. The stress-intensity factors calculated from the 
COD method assuming plane strain vere about 10 percent higher than those 
obtained assuming plane stress. Results from the nodal-force method fell 
generally between the two COD resvilts. Because neither plane-strain nor 
plane-stress conditions are fully realized in a cc.rner-cracked plate, either 

assumption is expected to lead to some inaccuracy in stress-intensity factors. 

/:< 

Quarter-Elliptical Corner Cracks /> 

Figures 10 and 11 show the normalized stress-intensity factors for 
q.uarter-elliptical corner cracks (a/c ~ 0,2 and 2, respectively) as a 
function of the parametric angle, (|), and the crack-depth-to-plate- 'j 
thickness ratio, a/t, for remote tension, remote bending, or wedge 
loading in the hole. 

Figure 10 (a/c = 0.2) shows that for remote tension, the larger a/t 
ratios gave higher normalized stress-intensity factors, but for remote 
bending or wedge loading the trend was generally reversed. For remote ^ 
bending, the maximum normalized stress-intensity factor occurred near 
(() = it/ 2 for a/t = 0.2 and at <() = 0 for a/t = 0.8. The wedge-loaded hole 
results (Fig. 10c) showed a steep gradient in stress-intensity factor near the 
hdle surface. 

Figure 11 (a/o = 2) shows that for remote tension and bending the 
smaller a/t ratios gave higher normrlized stress -intensity factors. For 
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remote t<?4ision and wedge loading, the stress-intensity factors varied 
little (about 10 to 15 percent) along the crack front compared to those 
for the other crack configurations (a/c » 0.2 and l). For remote bending, 
the normalized stress-intensity factor was largest^at <|» = 0 (front surface). 
Again, the negative stir ss-intensity factors shown in Figure lib for a/t = 
0.8 are meaningful only in the presence of sufficient tensile loading to 


prevent crack-surface contact. 

The stress -intensity factors (Mode I) for q.Uarter-elliptical corner 
cracks subjected to simulated pin loading, as shown in Figure 5, can be 
calculated by substituting into equation the results shown in Figures 10a 

and 10c for a/c = 0.2, and Figures 11a and 11c for a/c = 2. Figure 12 

' !i 

shows the normalized stress-intePolty factors for quarter-elliptical corner 
' Cracks (a/c = 0.2, 1.0, and 2.0)! as a function of the parametric angle, tj), 
and the crack-deptb-to-plate-thil'kness ratio, a/t, for simulated pin load- 
ing. The graphs are similar to those shown for the wedge-loaded case, 
differing mainly in magnitude. This similarity is due to the small R/h 
ratios considered here, which caused the contributions from the remote 
tension case to be small compared with those from the wedge-loaded case. 

As in the wedge-loaded case, for quarter-elliptic corner cracks with 
a/c = 0.2 the stress-intensity factors rise sharply with increasing <|) near 
the intersection of the crack with the hole surface. For quarter-circular 
corner cracks, the maximum stress-intensity factor occurs near (j> = tt/2. 

For quarter-elliptic cracks with a/c = 2, the stress-intensity factors 
are relatively constant along the crack front. 


CONCLUDING REMARKS 


A three-dimensional finite -element elastic stress analysis vas used 
to calculate Mode I stress-intensity factor variations along the crack 
front for two symmetrical q^uarter-elliptical corner cracks at the edge 
of a hole in a finite -thickness plate. The plates were subjected to 
remote uniform tension, remote bending, wedge loading in the hole or 

V' 

simulated pin loading in the hole. The ratio of crack depth to plate 
thickness ranged from 0.2 to 0.8, and the ratio of crack depth to crack 
length ranged from 0.2 to 2, while the ratio of hole radius to plate 
thickness was held at 0.5. Three-dimensional singularity elements in 
the shape of pentahedrons were used at the crack front-, elsewhere, eif^t- 
noded hexahedrons were used. A nodal-force method which requires no prior 
assumption of either plane stress or plane strain was used to evaluate 
the stress-intensity factors, 

A convergence study on stress-intensity factors for a deep quarter- 
elliptical corner crack in a finite-thickness plate showed satisfactory 
convergence. About 9300 degrees of freedom were used in calculating 
stress-intensity factors for all corner-crack configurations. ; 

For quarter-circular corner cracks subjected to tension or wedge 

loading in the hole, the maximum stress -intensity factor occurred near 

/' . .. 

the intersection of the crack with the hole surface, whereas for bending 
it occurred at the intersection of the crack with the front surface. 

For quarter-elliptical corner cracks (crack-depth-to-crack-length 
ratio of 0.2) subjected to tension or wedge loading in the hole, the 
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maximum stress-intenulty factor occurred near the intersection of the 
crack vith the hole surface for shallov cracks (crack-depth-to-plate- 
thickness less than or equal to 0.5) but occurred at the front surface 
for deep cracks' (crack-depth-to-plate-thickness ratio of 0.8). 

The stress-intensity factors for quarter-elliptical corner cracks 
(crack depth-to-crack length ratio of 2} were nearly constant (within 
about 10 percent) along the crack front in the cases of tension and wedge 
loading in the hole, though not in the case of bending, where the stress 
intensity factor was greatest at the intersection of r he crack with the 
front surface. 

For quarter-elliptical comer cracks subjected to simulated pin 
loading in the hole (balanced by uniform tension at one end) , the 
stress-intensity variations along the crack front were similar to those 
obtained for wedge loading in the hole, differing mainly in magnitude. 

The stress-intensity factors obtained herein should be useful in 
correlating fatigue-crack-growth rates and in calculating fracture toughnes 
for the corner-crack configurations considered. 
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TABLE I; Stress-Intensity Boundary-Correction Factors , F. , for 
Quarter-Elliptical Corner Cracks Subjected to Rlmote 
Tension, Remote Bending, and Wedge Loading in the Hole 

} _ 

(Fi = Kj/S^y7ra/Q ; b/(R+c)>5; h/b = 1.8; V= 0.3). 


TYPE 

OF 

LOADING 

a/c 

0.2 

jV- 

1 

2 


0.2 

0.5 

0.6/ 

0.2 

0.5 

0.8 

0.2 

0.5 

0.8 


0 

0.673 

0.878 

1 . 256 

1.832 

1.699 

1.808 

1.603 

1.353 

1.259 


0^125 

0.711 

0.934 

1.255 
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Fig. 5.“ Corner-crack configuration subjected to simulated 
pin loading. 
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Fig. 7.- Distribution of stress-intensity factors along crack front for 
quarter-circular corner cracks (a/c = 1; R/t = 0.5). 
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Fig. 8.- Comparison of stress-intensity factors for quarter- 
circular corner cracks subjected to remote tension 
(a/c = 1; R/t = 0.5). 
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Comparison of stress-intensity factors determined by 
displacement method (plane stress and plane strain) 
and nodal- force method (a/c = 1; a/t = 0.2; R/t = 0. 
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Fig. 10.- Distribution of stress-intensity factors along crack front for 
quarter-elliptic corner cracks (a/c = 0.2; R/t = 0.5). 
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Fig. 1’ - Distribution of stress-intensity factors along crack front for 
quarter-elliptic corner cracks (a/c = 2; R/t = 0.5). 









